If neutrinos are Dirac particles, then there wiU exist right-handed chirahty 
states of neutrinos un and left-handed chirality states of antineutrinos Vl which 
in the Standard Model Lagrangian do not have couplings to the weak interaction 
gauge bosons W and Z. The physical neutrino states are the helicity eigenstates iy± 
and V± , and because helicity and chirality eigenstates do not coincide (only for 
ultra-relativistic neutrinos ~ and ~ z^i?), all the neutrino and antineutrino 

(N ■ 

0> ' helicity states couple to the weak gauge bosons. Although the rate for W and 

ON ; 

Z mediated scattering between ultrarelativistic fermions involving one (i^-) is 
Q I suppressed, depending on the process at hand, by a factor (10~^ — l)(m^/2s), 
where s is the centre-of-mass energy squared, as compared with the corresponding 
^ I reaction involving u- (z^+), the rate might be sufficiently large to thermalize these 
'wrong-helicity' states in the early Universe. 

If neutrinos with 'wrong' helicities should still be present with an energy 
density comparable to at the time when neutron-to-proton ratio froze out (at 
T ~ 0.7 MeV), they would have speeded up the expansion rate of the Universe 
(which depends on the total energy density), with the consequence that the 
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limits. These imply that the effective number of extra relativistic two-component 
neutrino species ^A'',^ is much less than one, with SN^, < 0.3 perhaps the best 
limit [1]. Thus the production rate for z^_|_, and hence their mass [2], must 
be small enough so that they decouple already before the QCD phase transition 
which takes place at temperatures somewhere between 100 and 400 MeV. In that 
case they will not participate in the entropy transfer from quark-gluon plasma to 
particles in equilibrium, and consequently their number and energy densities will 
be diluted below levels that are acceptable for primordial nucleosynthesis. The 
production rate was first estimated by Fuller and Malaney [3], who argued that a 
Dirac neutrino with a lifetime exceeding the nucleosynthesis time scale (t ~ 1 s) 
should have a mass less than about 300 keV. 

This argument applies to neutrinos with a mass much less than 1 MeV. A 
heavy neutrino with a mass in the MeV region would have a more pronounced 
effect on nucleosynthesis than a light neutrino, because during the synthesis of the 
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light elements the energy density of the 'right-helicity' states of a heavy neutrino 
would be comparable to or higher than that of a massless neutrino. This is 
because at that epoch these states have already decoupled (T^g" ~ few MeV). This 
has been shown [4] to lead to an excluded region 0.5 MeV < m,^^ < 30 MeV for 
the tau neutrino mass, provided Tj^^ > 10^ s. (If 1 s < Tjy^ < 10^ s, the upper 
bound is somewhat weakened.) 

These considerations are relevant not only for the tau neutrino mass, the 
laboratory limit on which is m.,y^ < 31 MeV [5], but possibly also for the muon 
neutrino. This is because the exprimental mass limit has recently been revised 
upward [6], with the current limit being about m,y^ < 500 keV. 

The original Fuller and Malaney limit was based on an approximate estimate for 
the rate of production. In this Letter we shall present a careful re-evaluation of 
this limit by computing all the relevant cross sections and decay rates exactly. As 
we shall show by explicit calculation, the mass limit on Uj- (u^) becomes larger as 
much as by a factor of four (two) as compared with the Fuller and Malaney result. 

Before the QCD phase transition but below, say, T ~ 0.5 GeV, the fermions 
present in the Universe at significant number densities were the leptons and u, 
d, s and c quarks. All the 2 — > 2 scattering processes involving them, with no 
'wrong-helicity' neutrinos in the initial state and with at least one vf^ or z/^^ in 
the final state, are listed in Table 1. (There are altogether 47 separate reactions 
for each of the neutrinos which need to be taken into consideration). The first 
constraint is imposed on because each wrong helicity neutrino in the initial (final) 
state introduces an additional small factor m^/|pp (m^/|p'p) to the cross section. 
Here |p| and |p'| are the absolute values of the centre-of-mass momenta of the 
incoming and outgoing particles, respectively. Hence processes with more than one 
'wrong-helicity' neutrino can be ignored as compared to prosesses with only one 
'wrong-helicity' neutrino. 

We have computed the cross sections for these scatterings in the limit of low 
momentum transfer Iq^I ^ ^h- The contributions to the cross sections 

arising from the exchange of gauge bosons are given by 




(1) 
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where we have summed over all helicity and colour states of all charged leptons 
and quarks. The dimensionless functions F| ' are listed in Table 1. In Eq. 
(1) the coefficient of F^'^~'^'^^ is of the order of an ordinary weak interaction 
cross section for fermion-fermion scattering, so that all suppression due to the 
'wrong-helicity' neutrino in the final state is included in this function. For a 
massless neutrino ^j^^^^^^ = 0, and F^"^^^^"* increases monotonically as a function 
of its arguments. Some of the processes, such as u'^u'i u^u^ and z/^j^:L — > u^u'!^ 
{k = fj,,T ^ j = e, fj,, r) require two mass insertions and are therefore suppressed. 
Similar suppression holds in the relativistic limit also for the charged-current 
processes and l^~u'^ v^iF' {u^ = u,c; = d,s), despite the 

fact that there is only one 'wrong-helicity' neutrino. 

Having exact formulas one can estimate the cross sections for 'helicity-flip' 
scattering in the relativistic limit. It turns out that the popular approximation 
a+ ^ G']pEl{mjy /2EjyY , where the overbar indicates the averaging over the helicity 
and colour states of the initial charged leptons and quarks, overestimates the cross 
section for all the processes typically by more than an order of magnitude. 

Neutral current scattering between fermions can also be mediated by the Higgs 
boson H. Although the Yukawa vertices provide direct 'helicity-flip' interactions 
for fermions, the fermion couplings to Higgs are weaker than gauge couplings 
roughly by a factor ruf/Mw- In the low energy limit additional suppression would 
arise from the propagators if the mass of the Higgs is larger than Mw and Mz- 
We have calculated the contributions from the Higgs boson exchange and from the 
interference between gauge boson(s) and the Higgs boson, and found that these 
are smaller at least by a factor |p'p/M^ than gauge boson contributions to the 
cross section. In what follows we shall neglect the Higgs boson contribution. 

Besides the scattering processes, and Uj- can be created also in three-body 
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decays. There exist six relevant decay channels producing vi^: 
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and another six decay channels producing v 
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v'l^du, 
u^dc, 



zv^s c, 



The decay rate for particles, e.g. —>■ v'^e Ug, is in the rest frame given by 
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(mi - v^|p4p + m|)(miA/|p4p +m^ -mi|p4| - m|) 
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(V|P4p - |p4|) 

while the decay rate for antiparticles, e.g. 17^//"'" i^^, is given by 

|P4p(^P4pTmf - |p4|) Ai/2(m?4, ml, m|) 



.(1^234) 



{2A(mi4,m2,m3) + [mi4 + [ml + ml)ml^ - 2{ml - mlf]], 



m 



14 



with = mf —2m i-\/jp4P"TmJ +777.4 A(a, b, c) = a^ + 6^ + c^ — 2a6 — 26c — 2ca. 
Here the subscript 4 refers to or z/^ and the iiiaxinial nionicntum of the 'wrong- 
hehcity' neutrino is |p4|maa; = ^/[mf — (m2 + ma)^ + ml]^ — 4mi^m42/2mi. The 
expressions for particle and antiparticle decay rates are different because they do 
not correspond to CP-conjugated processes: in both cases the final state involves 
the 'wrong-helicity' neutrino (but not anti-neutrino). In these rates we have 
summed over all helicity states of all particles except the u^^'^^ under consideration. 
If quarks are included in the decay process these decay rates must be multiplied 
by a factor SV^^ which accounts for the three colour states of quarks and for 
quark mixing. 

The thermally averaged scattering rate reads [7] 

rf = ^ ^]:J ^.^JiE.mfiE.m.T~-'^i,un), (5) 

where n^^ is the equilibrium number density of z^+'s, f{Ei/T) are the Fermi-Dirac 
distributions of the incoming particles, and j{pi,P2) = a/ (pi • ^2)^ — mfm|/E'i£'2 
is a flux-related factor. In (5) we have neglected the final state Pauli blocking, 
which is an about 10% effect [7]. 

Thermally averaged decay rate is simply given by 

where the factor mi/Ei arises from the Lorentz boost of the decay rate. 

We have estimated the total thermally averaged z/_|_ production rate 
r+ = r^ + T'^ numerically, and the result is displayed in Fig. 1. The difference 
between the production rates for and u]^ is due to the differences in the phase 
spaces (initial and/or final) of the corresponding processes. While fj, is relativistic 
in the range of temperatures 100 — 400 MeV, r is nonrelativistic and therefore 
its number density, being Boltzmann-suppressed, is very sensitive to changes in 
temperature. This can be illustrated by the fact that while at T fa 100 MeV purely 
charged processes do not contribute to the production, their contribution is 
about 40% at T 250 MeV and about 75% at T ^ 400 MeV. In contrast, for 



this figure rises from 55% (T a; 100 MeV) to 70% (T fa 400 MeV). One can see 
that elastic scattering and annihilation processes, dominated by i/^I7^ — > i/^l7^, 
are not so important as one would naively expect. Comparing with the Fuller and 
Malaney result, which is also shown in Fig. 1, we find that at T Ri 100 MeV the 
actual rate for {ly'^) production is 18 (7) times smaller. 

We now require that the production of z^_|_'s ceases latest at the onset of the 
QCD phase transition, or that 

r+(m„TQCD) <iy(TQCD), (7) 
where the Hubble parameter H is given by 



Here p(mj, T) is the total energy density, including also the particles that could 
be non-relativistic at the time of QCD phase transition, such as r, c, s and 
fi. Here we differ from the treatment of Fuller and Malaney, who did not 
account for the non-relativistic degrees of freedom in their estimate. Thus 
geff counts all the degrees of freedom, and we have tabulated it in Table 2, 
assuming that the non-relativistic species stay in equilibrium with the appropriate 
Boltzmann suppressed equilibrium densities. We have used the current quark 
masses m^, 5.6 ± 1.1 MeV, rud ^ 9.9 ± 1.1 MeV, 199 ± 33 MeV, and 

rric ~ 1.35 ±0.05 GeV [8]. About the half of the difference between the minimal 
and maximal values of Qeff is due to the uncertainty in quark masses. Another 
half comes from the uncertainty as to whether both or only one of the neutrinos 
u'^ and is in equilibrium. 

As the production rate for the 'wrong-helicity' neutrino increases together 
with its mass (at T » m,^ r+ oc ml), and for the temperatures and masses 
under consideration H does not depend on rui,, there exists an upper limit on the 
neutrino mass at which the inequality (7) can be satisfied. These mass limits for Vt- 
and as functions of the QCD phase transition temperature are displayed in Fig. 
2. The value of the upper limit of the neutrino mass is almost insensitive to the 
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variations of the quark masses within ranges given above, despite the fact that the 
number and energy densities of s and c quarks change considerably. This happens 
because the changes in the production rate and the expansion rate compensate 
for each other. From Fig. 2 one can readily see that if Tqcd — (100)200 MeV, 
the mass limits are < (1180)740 keV and m,.^ < (720)480 keV. According 
to ref. [4], however, such large masses are already in conflict with primordial 
nucleosynthesis. Thus we may conclude that the equilibration of 'wrong' Dirac 
neutrino helicity states does not yield any additional new limit. Note that because 
the dilution provided by the entropy production at the QCD phase transition is 
more than sufficient for nucleosynthesis, v+s could well decouple some time during 
the actual phase transition, lowering the decoupling temperature and increasing 
the cosmological neutrino mass limit. Assuming decoupling at the onset of the 
phase transition will therefore slightly underestimate the actual upper limit. 
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{l/12)x^k{VTT^+lf 
(l/6)x,.[2(4 T A? + 3(4' + 34^)0;,,] 
(l/6)a;,fe[2(4 c\f + 3(4 + Sc^^)^,^.] 
^.^9 ^ ^/^^ (l/2)x,.[2(4 - c\f + 3(cf? + 3c?)x,] 

z^^g ^ z.^g (l/2)x.. [2(4 + c\f + 3(4 + 'icf)x,\ 

^ z/^z7^ (l/12)y,. (v/TT^ + 1)2 

^;c-^fc+ ^^fc^fc (l/3)y,.(42 + cf )(2 + 3x^0 
^ ^,^27^ (l/3)y,.(4 + 42)(2 + 3xi,) 
^ y,, (4 + ) (2 + 3x,) 

ik-^3_ ^ ^fep- 2(v/rn^ - i)(v/rH^+ i)(v/rq^ + i){^/tt^ - 1) 

^ (2/3)(^/rT^- i)(A/rH^+ i)(3^/rq^v^rTy^- 1) 

^ j,;fc77^^ (2/3)(vi + 2/.fc -i)(v/T+y;7 + i)(3yrq^^/rq^ - 1) 

;fe-^,m^^fc^n i2v^,(v/i + y.'c -i)(vi + yd" -i)(v/rT^vi + ^K"^ + 1) 

4v;^„(Vi + y^^ - i)(3v/rT^Vi + ^^d- Vi + Vu^ - i) 
^ 4v;^^(^m/;;^ - i)(3vi + xum^i + VT+j^ - 1) 



Table 1. List of 2 — > 2 scattering cross sections for production. Fj^ ~^ 
is defined by Eq. (1), and the notations are: k — n, r ^ j — e, //, r; = 
u, c; (T^ = d, s; Xf = m^/|pp and yf = m^/|p'p; 4 = + 1) = + 1) 
Cy = T^^ — 2Q^sui'9w and c;^ = T^-'^ are the vector and axial vector couplings; 
Vmn is an element of the C-K-M matrix. We have used the values sin' 9w = 0.2325, 
Vud = 0.9753, Vus = 0.221, = 0.221 and Vcs = 0.9743 [8]. The present 
experimental uncertanties in these values do not affect our upper limits on the 
neutrino masses. 
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rTi/-\ jr^\ T\ 

1 (JVLeVj 


3eff 


nia,x 
3eff 


100 


59.4 


62.8 


150 


61.6 


64.3 


200 


62.8 


65.3 


250 


64.0 


66.5 


300 


65.3 


67.8 


350 


66.7 


69.1 


400 


67.9 


70.3 



Table 2. The effective number of degrees of freedom as a function of temperature. 
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Figure Captions 



Fig. 1. The total thermally averaged production rate for ly^ and i^^. For 
comparison, we show also the Fuller and Malaney rate (dashed curve). 

Fig. 2. Neutrino mass limits corresponding to r+ < H. The allowed region is 
below the curves. 
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Abstract 

We consider the equihbration of the 'wrong-hehcity' Dirac neutrino states and 
T'_ in the early Universe via weak interactions and calculate carefully the thermally 
averaged production rate, taking into account all the relevant scattering and decay 
processes. Requiring that the production rate is less than the Hubble parameter 
at the onset of QCD phase transition so that the nucleosynthesis predictions are 
not contradicted, we find for Tqcb — 200 MeV the upper limits m^^ < 740 keV 
and < 480 keV. 
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